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We present tunneling and rattling motions of an off-center guest atom in cage referring to a 
clathrate crystal LasPcboGeg. The elastic constant C44 of LaaPcboGee shows a Debye-type disper- 
sion around 20 K obeying a relaxation time r = roexp(_E / ksT) with an attempt time to = 2.0 x 10~ 12 
sec and an activation energy E = 197 K. At low temperatures below 3 K down to 20 mK, the C44 
shows a softening of C* 44 = C° 4 (T - T%)/(T - 9) with ig = -337.970 mK and = -338.044 
mK. These facts are attributed to two different types of the off-center motions with Ps symmetry 
in 4a-site cage of La3Pd2oGe6, a thermally activated rattling motion over the potential hill and a 
tunneling motion through the potential hill at low temperatures. 



PACS numbers: 63.20.Pw,62.20.Dc,82.75.-z 

I. INTRODUCTION 

It has recently been found that filled skutteruditcs 
RM 4 Sbi2 (R: La or Ce, M: Fe or Co) [1], clathrate semi- 
conductor compounds SrgGaigGeao [2] and EusGai6Ge3o 
[3] show considerable reduction of the thermal conduc- 
tivity. The rattling motion of an off-center guest atom 
over a potential hill in oversized cage in these clathrate 
compounds brings about reduction in thermal phonon 
transport by a resonant phonon scattering, which is fa- 
vorable for thermoelectric devices with profitable figure 
of merit [4]. Furthermore, the low-temperature behavior 
of thermal conductivity k(T) ~ T 2 and ultrasonic at- 
tenuation a(T) ~ T 3 in Sr 8 Gai6Ge3o [2] reminds us of 
the two-level system due to tunneling motions in glasses 
with structural disorder [5-8]. It should be noted that 
the off-center guest atoms of clathrate compounds is lo- 
cated in cages forming an ideal periodic lattice. The off- 
center tunneling in clathrate crystals, therefore, may be 
regarded as a new quantum degrees of freedom associ- 
ated with a local Einstein phonon with a large density of 
state. 

Quite recently, our group has found the frequency 
dependence in the elastic constant C44 of Ce3Pd2oGe 6 
showing a ferroquadrupole ordering [9,10] and in 
(Cn — Cxi)/2 of a filled skutterudite PrOs4Sbi2 show- 
ing a heavy fermion superconductivity [11]. The ultra- 
sonic dispersion in these compounds shows the fact that 
a relaxation time for the rattling motion obeys thermal 
activation-type temperature dependence. We have suc- 
cessfully projected out the off-center mode in cage, which 
couples to appropriate ultrasonic mode. The phenomena 
associated with the tunneling motion are expected at low 
temperatures, where the thermally activated rattling mo- 
tion dies out completely. However, the ferroquadrupole 
ordering at Tq = 1.25 K in Ce3Pd2oGe6 and the su- 
perconductivity at Tc = 1-85 K in PrOs4Sbi2 smear out 
the characteristic low-temperature properties due to tun- 
neling motions. In the present paper, we have made 



ultrasonic measurements on La3Pd2oGeg with a cubic 
Cr23C6-type structure in order to examine the tunnel- 
ing and rattling motions in the system being absent of 
the 4f-electron. 



II. EXPERIMENT 

An image furnace equipped with two ellipsoidal mir- 
rors was employed to grow a single crystal of La3Pd2oGe 6 
with residual resistivity po — 2.7/if2cm. The plane par- 
allel surfaces with (100) or (110) orientation were deter- 
mined by x-ray photograph. The sound waves were gen- 
erated and detected by the piezoelectric LiNbC>3 trans- 
ducers bonded on the surfaces of specimen. The change 
of sound velocity v was detected by a phase comparator 
based on mixer technology. The elastic constant C = pv 2 
was estimated by the mass density p = 10.179 g/cm 3 of 
La3Pd2oGe6 with lattice parameter a = 12.482 A [11]. A 
3 He-evaporation refrigerator down to 450 mK and a 3 He- 
4 He dilution refrigerator down to 20 mK were employed 
for low-temperature ultrasonic measurements. Observa- 
tion of the de Haas oscillation in elastic constants [12] 
ensures good quality of the present single crystal. 

III. RESULTS AND DISCUSSIONS 

We show temperature dependence of elastic constants 
C u , (Cn - Ci 2 )/2 and C 44 of La 3 Pd2 Ge 6 in Fig. 1. 
Here the longitudinal sound wave with frequencies of 
70 MHz was employed for measurements of Cn and 
transverse one with 40 MHz for (Cn — Cx-x)j1 and C44. 
With lowering temperature we have found a step-like in- 
crease around 20 K in C44, while other modes of Cn and 
(Cn — Ci2)/2 exhibit monotonous increase only. 

The clastic constant C44 of transverse acoustic mode in 
La3Pd2oGe6 in Fig. 2(a) exhibits a marked dependence 
in frequencies for 9 MHz, 40 MHz and 150 MHz. This 
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FIG. 1: Temperature dependence of elastic constants Cn, 
(Cii — Ci2)/2 and C44 of La3Pd2oGe6- A step-like anomaly 
in C44 around 20 K indicates the ultrasonic dispersion. 
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result indicates the ultrasonic dispersion in the C44 mode 
as similar as that in isomorphous compound Ce3Pd2oGe6 
[9]. These results show the fact that R 3 Pd2oGe6 (R = 
La, Ce) exhibit commonly a thermally activated rattling 
motion of the off-center rare-earth ion. As was noted in 
Fig. 1, the ultrasonic dispersion has been found only in 
C44 , but no sign of dispersion in {C\\ — C\2)/2 and Cn. 
This result immediately means that the off-center motion 
has the F 5 symmetry. Consequently, the charge fluctua- 
tion due to the off-center motion couples to the acoustic 
transverse C44 mode. On the other hand, a filled skut- 
terudite PrOs4Sbi2 of heavy fermion superconductor ex- 
hibits ultrasonic dispersion in (Cn — Ci2)/2 indicating 
the T 2 3 symmetry of the off-center mode in cage [11]. 
It is worthwhile to emphasize that the symmetry trans- 
verse ultrasonic mode is a useful probe to determine the 
symmetry of the off-center motion in cage. 

The frequency dependence of the elastic constant in 
Fig. 2(a) is described by the Debye-type dispersion as 

r (, A r ^44(00) -C 44 (0) m 
644(0;) = G 44 (oc) — — — . 1 

1 + LO Z T Z 

The calculation of inset of Fig. 2(a) reproduces well 
the experimental results in Fig. 2(a). Here lo is the angu- 
lar frequency of the sound wave. The relaxation time r 
of the rattling motion in Fig. 2(b) revealed the ther- 
mally activated-type temperature dependence. 6*44 (0) 
and 644(0x0) are low- and high-frequency limits, respec- 
tively. Arrows in Fig. 2(a) indicate the temperatures 
being satisfied with the resonant condition, where the 
relaxation time r coincides with the sound frequency to 
as lot — 1. The ultrasonic attenuation of the transverse 



FIG. 2: (a) Frequency dependence of the elastic constant C44 
measured by the transverse ultrasonic waves with frequencies 
of 9 MHz, 40 MHz and 150 MHz. Arrows in (a) indicate 
the temperatures satisfying with the resonance condition of 
lot — 1. Inset in (a) represents a fit by the Debye-type dis- 
persion of Eq. (1) for frequencies of 9 MHz, 40 MHz, and 
150 MHz. (b) shows the Arrhenius plot of relaxation time 
t — Tocxp(E / ksT) in La3Pd2oGe6 together with the results 
of Ce 3 Pd2oGe 6 and PrOs 4 Sbi 2 [9,11]. 



644 mode, that is not presented here, shows a peak at the 
temperatures satisfied with lot — 1 . This resonance scat- 
tering of the ultrasonic C44 wave by the rattling motion 
in cage may bring about a reduction of thermal conduc- 
tivity by the phonon transport. 

It should be noted that the specified symmetry 
of the transverse C44 wave may give rise to anisotropic 
phonon transport depending on the crystallographic ori- 
entation. The longitudinal Cn mode as well as the 
transverse C44 mode participates in the thermal phonon 
transport along the [100] direction. The later C44 
mode is scattered by the rattling motin, while the for- 
mer Cn mode propagate transparently without any 
scattering by the rattling. In the case of thermal 
phonon transport along the [111] direction, the longi- 
tudinal (Cn + 2C12 + 4C44)/3 as well as the transverse 
(Cn — 2C12 + Ct4)/3 are scattered by the rattling, be- 
cause both modes contain the pure C44 mode in part 
[13]. These anisotropic properties of acoustic phonon 
may bring about the anisotropy of the thermal conduc- 
tivity by phonons. The investigation of the ratting by 
ultrasonic measurements would be profitable way to de- 
velop the thermoelectric device with high figure of merit. 
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As one can see in an Arrhenius plot of Fig. 
2(b), the relaxation time r of La 3 Pd 2 oGe 6 obeys 
the thermal activation-type temperature depen- 
dence r = ToCxp(E/kBT) with an attempt time 
To = 2.0 x 10~ 12 sec and an activation energy E = 197 
K. The relaxation times r in Ce3Pd2oGe6 with 
t = 3.1 x 1CT 11 sec, E = 70 K [9] and in PrOs 4 Sbi 2 
with t = 8.8 x 10" 11 sec, E = 168 K [11] are also 
presented in Fig. 2(b) for comparison. These relaxations 
originate from the thermally assisted rattling motion 
of the off-center rare-earth ion over a potential hill 
with height E = 70 ~ 200 K. It is of importance to 
contrast with the charge fluctuation time of 4f electron 
in inhomogeneous mixed valence compound Siri3Te4 
showing the activation-type temperature dependence 
with T = 2.5 x 10~ 13 sec and E = 1600 K [14]. The 
relatively slow relaxation time To = 0.2 ~ 8.8 x 10~ n 
sec as well as relatively small activation energy 
E = 70 ~ 200 K in La 3 Pd 20 Ge 6 , Ce 3 Pd 20 Ge 6 and 
PrOs4Sbi2 indicate the rattling motions of the off-center 
rare-earth ions with heavy masses. The motion of 
rare-earth atom described by a harmonic oscillation of 
C(z) — (l/7r^o) 1 ^ 2 exp(— z 2 /2zq) in cage spreads in space 
with a mean square displacement zq = (l/2ir)(hTo/M) 1 / 2 
[15]. The attempt time to = 2.0 x 10~ 12 sec and the 
mass M = 139m p of La atom in La3Pd2oGe6, where 
m p is proton mass, lead to a mean square displacement 
zo = 0.12 A being comparable to zq = 0.48 A of 
Ce 3 Pd 2 oGe 6 , z = 0.79 A of PrOs 4 Sb 12 [9,11]. Further- 
more neutron scattering on EusGai6Ge3o shows zq = 0.8 
A [16]. 

In order to examine low-temperature properties due to 
tunneling motions of La3Pd 2 oGe6, we have made ultra- 
sonic measurements down to 20 mK. As shown in Fig. 
3, we have observed a remarkable softening of C44 with 
relative change 200 ppm below 3 K down to 20 mK by 
use of the dilution refrigerator. The (Cn — Ci 2 )/2 mode 
in inset of Fig. 3 shows a softening of 20 ppm, which is 
much smaller than the 200 ppm softening of C44. The 
characteristic softening in C44 of Fig. 3 in particular 
promises the T§ symmetry of the off-center tunneling 
motion in cage. The tunneling frequency is estimated 
about 60 GHz corresponding to the onset temperature 
3 K in the elastic softening of Fig. 3. The present ul- 
trasonic measurements with frequencies 10 ~ 200 MHz 
are able to observe the static strain susceptibility for the 
off-center tunneling mode. Magnetic fields up to 8 T do 
not make any appreciable change in the low-temperature 
softening of C44. This result, which is not presented here, 
ensures that the softening of C44 is intrinsic effect being 
independent of magnetic impurity effects. 

In the present clathrate compound La3Pd 2 oGe6 with 
space group symmetry Fm3m, there are two different 
types of cages containing La ion. The unit cell of 
La3Pd 2 oGe6 with a lattice parameter a — 12.482 A con- 
sists of four molecular units with 116 atoms. The 4a- 
site cage with Lai atom consisting of twelve Pd-atoms 
with distances G?Lai-Pd2 = 3.084 A and six Ge-atoms 



t 1 r 




12 3 4 

T(K) 



FIG. 3: The low-temperature softening in elastic constant 
C44 of La3Pd2oGe6 below 3 K down to 20 mK. Solid line is 
a fit by C44 = C2JT - Tg) / \T - 9) with the parameters in 
the text. Inset presents low-temperature behavior of C44 and 
(Cii - Ci»)/2. 



with o?Lai-Gc = 3.350 A forms a face center cubic lattice. 
The 8c-site cage with La2 atom consisting sixteen Pd 
atoms with G?La2-Pdi = 2.884 A and rfLa2-Pd2 = 3.392 
A forms a simple cubic lattice [10]. The trivalent La 
ion with radii 1.88 A in former oversized 4a-site cage 
may favor the rattling motion over off-center positions. 
Actually neutron scatterings on R,3Pd2oGe6 (R = Ce, 
Pr, Nd) observed well-defined crystalline electric field ex- 
citations of 8c site, while they detected obscure peak 
only for 4a site [17]. These facts are consistent with 
the off-center motion of the rare-earth ion in the 4a site 
cage of R.3Pd2oGe6. Because the absence of Pd- or Ge- 
atom on the three- fold [111] axis in the 4a-site cage, 
the potential for Lai ion may possess off-center minima 
at ri = (a, a, a), r-2 = {—a,— a, a), = (—a, a,— a), 
r 4 = (a, -a, -a), r 5 = (a, a, -a), r 6 = (-a, -a, -a), 
rj = (a, —a, a), r§ = (—a, a, a) with a ~ zo/2 = 0.06 A. 
Here r^(i = 1,2,..., 8) denotes distance from the center 
of the 4a-site cage. 

It is useful to project out the irreducible representation 
for the off-center modes consisting of the atomic densities 
Pi = p{ r i) a t the off-center positions r» in 4a-site cage 
with the site symmetry Oh [9, 18, 19]. The off-center 
motion over the eight off-center positions is reduced into 
direct sum T\ © T 2 (BT4 © r 5 . Because the C44 mode 
with r$ symmetry in La3Pd2oGeg shows the ultrasonic 
dispersion around 20 K in Fig. 2(a) and the appreciable 
low-temperature softening below 3 K in Fig. 3, the 
off-center mode prs.yz = Pi-p2-P3+Pi-pa+P6~Pi+Ps, 

PTS.zx = Pi - P2 + P3 - Pi - P5 + P6 + Pi - PS, PT$,xy = 

Pi + P2 - P3 - Pi + P5 + P& - Pi - Ps is relevant for the 
ground state in the present system. The T§ mode consists 
of anisotropic fractional atomic state over eight off-center 
positions rj(i = 1,2,..., 8). The prs.xy mode of Fig. 4, 
for instance, has fraction 1/4 at r\,r-2,r§, r§ positions and 
null at T3, T4, T7, r%. Other modes with T\, T2 or T4 sym- 
metry probably correspond to exited states and hence 



they are disregarded here. The Schrodinger equation for 
the off-center atom in multi potential wells is a faithful 
way to analyze the off-center tunneling state [20,21]. The 
present group theoretical analysis is a proper way instead 
of rigorous treatment. 

Because the La ion is trivalent, the charge density as- 
sociated with the r$ off-center mode in the 4a-site cage 
is described as Ap = Q yz prs, y z + QzxPrb.zx + QxyPrs.xy 
[18,19]. Here Q yz ,Qzx and Q xy denote the normal coor- 
dinate describing the T 5 off-center mode in Fig. 4. It is 
notable that the hexadecapole as well as quadrupole are 
relevant for the charge distribution of the off-center 
mode. As the temperature is lowered, the phase transi- 
tion Ap ^ due to the condensation of the symmetry 
breaking T5 mode would be expected in principle. We 
refer Yb4As3 to the phase transition due to charge or- 
dering [18]. The present system, however, no phase tran- 
sition has been observed down to 20 mK. Consequently, 
the tunneling motion of the off-center La-ion through the 
potential hill in keeping the Oh site symmetry is relevant 
even at low temperatures. 

The interaction of the T§ off-center tunneling mode in 
cage to the elastic strain e yz , e ZXl e xy of the C44 mode is 
described as [22-24] 

Hqs = ~5r 5 ^{Q y z(i)e y z + Qzx{i)£zx 

i 

+Q xy (i)e xy }. (2) 

Here grs is a coupling constant and means summa- 
tion over cages at site i. The inter-cage coupling of the 
T§ off-center tunneling mode accompanying the charge 
fluctuation is introduced as 

H QQ = -9r 5 ^2{(Qyz)Q yz (i) + {Qzx)Qzx{i) 

i 

+ {Qxy)Qxy(i)}- (3) 

Here g' T5 is a coupling constant and (Qij){i,j — x, y, z) 
denotes a mean field. The low-temperature elastic soft- 
ening in C44 below 3 K in Fig. 3 is described by 

O44 - G 44 - - j—. (4 ) 

Here N is the number of cages in unit volume. The 
strain susceptibility for the T 5 off-center mode in cage 
is written as gl 5 Xrb{T) = C Q /T = S 2 /T. Here is 
5 = (Cq) 1 / 2 is a deformation coupling energy for the 
r$ ground state. Consequently, one obtains C44 = 
CUT - T°)/(T - 0) with = (g' r5 /g 2 r5 )C Q and T» = 
NCq/C% 4 + 0. The solid line in Fig. 3 is a fit by the 
characteristic temperatures = —338.044 mK, TjS = 
-337.970 mK and a back ground = 3.33085 x 10 10 
J/m. The low-temperature softening of C44 propor- 
tional to reciprocal temperature in La 3 Pd 2 oGe6 is at- 
tributed to the Ts tunneling mode with the symmetry 
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FIG. 4: Schematic view of the Fs off-center mode pr 5 ,xy = 
Pi + P2 — pz — p4 + Ph + Po — P7 — ps possessing fractional states 
of Lal-ion over eight off-center potential minima along three 
fold [111] axis in 4a-site cage of La3Pd2oGe6- 



breaking character in the present system. This is strik- 
ingly different from minima in elastic constants of NaCl, 
which has an off-center tunneling motion with Ti singlet 
ground state of impurity OH ion [25,26]. 

The value of NCq/C^ =T% - 6 = 0.074 mK corre- 
sponds to the energy gain due to the coupling of the 
T5 off-center mode to the elastic strain as similar as 
the Jahn- Teller energy in d- or f-electron system [22- 
24]. Employing the number of the 4a site in unit vol- 
ume, N — 2.057 x 10 21 cm~ 3 , one obtains the defor- 
mation coupling energy S = 9.3 K corresponding to the 
energy shift for the T§ mode by unit external strain. It 
is remarkable that the deformation coupling of about 10 
K for the present tunneling mode is order of magnitude 
smaller than that of 100 K in 4f-electron system of rare- 
earth compounds and of 1000 K in 3d-electron system of 
transition metal compounds. This means that the charge 
fluctuation due to the tunneling motion is well screened 
by the cage. The negative value of = —338.044 mK 
indicates the antiferro-type inter-cage interaction of the 
T§ mode among the 4a-site cages with nearest neighbor 
distance 8.8 A. 

The tunneling motion in the two-level system of many 
structural glasses [5,27] provides the principal concept to 
elucidate the off-center tunneling motion of the clathratc 
crystals. The two-level system in ramdom potential of 
structural glasses losses its symmetry property. This 
is strikingly different from the specific symmetry of 
the off-center tunneling motion in periodically arrayed 
cages of the clathrate crystal La3Pd2oGeg. As was al- 
ready shown, the off-center tunneling in the present 
La3Pd2oGe6 has small antiferro- type inter-cage interac- 
tion = —338.044 mK. Therefore, the collective motion 
of the off-center tunneling in periodically arrayed cages is 
expected to show a fiat dispersion being mostly indepen- 
dent of the wave vector. The triply degenerate tun- 
neling state may be regarded as a local Einstein phonon 
with large density of state. 

It is expected that the off-center tunneling motion in 
cage of La3Pd2oGe6 interacts to the conduction electron 
consisting of the 4d orbit of Pd and 4p orbit of Ge of 
cage. This interaction of the local Einstein phonon to 
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the electrons may play a crucial role to bring about the 
small activation energy E = 197 K and the apprecia- 
ble mean square displacement z = 0.12 A. We noted 
many theoretical investigations showing that in a strong 
electron-phonon coupling regime, the renormalized po- 
tential for a guest atom may favor off-center positions 
and reduction of the excitation energy [15, 28-29]. The 
screening of the off-center tunneling mode, namely local 
Einstein phonon is mapped to the multi-channel Kondo 
model of the non-Kramers doublet in f-electron systems 

In the case of the heavy Fcrmion superconductor 
PrOs4Sbi2, the ultrasonic dispersion has been found in 
the (Cn — Ci2)/2 mode [11]. The doubly degenerate 
T23 tunneling state accompanying the charge fluctuation 
may couple with the quadrupole fluctuation of the IV 
T\ ' pseudo quartet CEF state as well as the conduc- 
tion band [31]. This interaction of the local phonon to 
the electrons is of crucial importance for the formation 
of the heavy Fermion state and its Cooper pair of the 
superconductivity. The tunneling state with doubly or 
triply degeneracy in cage of clathrate compounds may 
regarded as a new quantum degrees of freedom leading 
to exotic low-temperature properties concerning the sym- 
metry breaking phase transition, superconductivity, and 
the multi-channel Kondo effect. 



IV. CONCLUSION 

In conclusion, the present ultrasonic measurements on 
the clathrate compound La 3 Pd 2 oGe6 revealed the dis- 
persion in the transverse C44 mode around 20 K due 
to the thermally assisted rattling motion over the po- 
tential hill. We have also found the appreciable elastic 
softening of C44 below 3 K down to 20 mK attributed 
to the off-center tunneling motion through the poten- 
tial hill. The off-center motion in the present clathrate 
crystal La3Pd2oGe6 is an ideal local Einstein oscillator 
with triple degeneracy in periodically arrayed cages. The 
clarification of the exotic properties of the local Einstein 
phonon being strongly coupled to electrons in metallic 
clathrate compounds is an interesting issue of strongly 
correlated physics in future. 
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